Quantum ratchets in dissipative chaotic systems 
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Using the method of quantum trajectories we study a quantum chaotic dissipative ratchet ap- 
pearing for particles in a pulsed asymmetric potential in the presence of a dissipative environment. 
The system is characterized by directed transport emerging from a quantum strange attractor. This 
model exhibits, in the limit of small effective Planck constant, a transition from quantum to classical 
behavior, in agreement with the correspondence principle. We also discuss parameter values suitable 
for implementation of the quantum ratchet effect with cold atoms in optical lattices. 
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Cold atoms exposed to time-dependent standing waves 
of light provide an ideal test bed to explore the features 
of the quantum dynamics of nonlinear systems. They 
allowed the experimental investigation of several impor- 
tant physical phenomena, such as dynamical localization 
[1, 2], decoherence [3], quantum resonances [4], and chaos 
assisted tunneling [5, 6]. In particular, they made pos- 
sible to realize experimentally the quantum kicked rotor 
[1, 2, 3, 4], a paradigmatic model in the field of quantum 
chaos [7]. Furthermore, cold atoms in optical lattices al- 
lowed to demonstrate the ratchet phenomenon [8, 9, 10], 
that is characterized by directed transport in the absence 
of any net force. 

Directed transport [11, 12] is found in periodic sys- 
tems due to a broken spatial-temporal symmetry [13], 
for instance in presence of lattice asymmetry, unbiased 
periodic driving and dissipation. This phenomenon, also 
called ratchet effect, is of potential relevance in tech- 
nological applications such as rectifiers, pumps, parti- 
cle separation devices, molecular switches, and transis- 
tors. Moreover, it is of great interest for the under- 
standing of molecular motors in biology [14]. In previous 
works, attention has been focused on systems with ex- 
ternal noise, on chaotic dynamical systems with dissipa- 
tion and on purely Hamiltonian systems (with mixed or 
chaotic phase space) [15, 16, 17, 18, 19, 20, 21]. In par- 
ticular, the issue of current reversal has attracted much 
interest [22, 23, 24, 25]. It was found that current rever- 
sal can be originated by a bifurcation from a chaotic to a 
periodic regime, that is, by the transition from a strange 
attractor to a simple one [23] . 

The presence of strange attractors is a characteristic 
feature of classical dissipative chaotic systems [26]. In 
quantum mechanics, it was found that the fractal struc- 
ture of the classical strange attractor is smoothed on 
the scale of Planck's cell [27]. It is therefore interest- 
ing to investigate how this phenomenon affects quantum 
ratchets. Till present the effects of dissipation have been 



mainly analyzed for classical ratchets while for quantum 
motion the theoretical studies have been concentrated 
on the Hamiltonian case [17, 18, 19]. In this Letter, 
we study a quantum ratchet system which is in contact 
with a dissipative environment and whose classical coun- 
terpart is chaotic. The investigation is carried out for 
parameter values where a strange attractor is present. 
Using the quantum trajectories approach, we are able 
to simulate numerically the evolution of a fractal quan- 
tum ratchet emerging from a quantum strange attractor. 
Moreover, we show that the ratchet effect is robust in 
presence of noise. Finally, we consider parameter values 
which should allow the observation of quantum chaotic 
dissipative ratchets in cold atoms experiments. 

We study a particle moving in one dimension [x G 
(— oo, +oo)] in a periodic kicked asymmetric potential: 

+ oo 

V(x,r)=k cos(x) + ^ cos(2x + </>)] 5(r-mT), 

m— — OO 

(1) 

where T is the kicking period. The evolution of the sys- 
tem in one period is described by the map 



n = jn + k(sm(x) + asin(2a; + 0)), 
x — x + Tn, 



(2) 



where n is the momentum variable conjugated to x and 7 
is the dissipation parameter, describing a velocity propor- 
tional damping. We have < 7 < 1; the limiting cases 
7 = and 7=1 correspond to overdamping and Hamil- 
tonian evolution, respectively. Introducing the rescaled 
momentum variable p = Tn, one can see that classical 
dynamics depends on the parameter K = kT (not on k 
and T separately). We note that at a = the classi- 
cal model reduces to the Zaslavsky map [29]. However, 
due to the symmetry of the potential, the ratchet ef- 
fect is absent in this case (see discussion below). The 
quantum model is obtained by means of the usual quan- 
tization rules [7]: x — > x, n — > h = —i(d/dx) (we set 
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ft = 1). Since [x,p] — iT, the effective Planck constant 
is Ti e ff = T. The classical limit corresponds to h e s - * 0, 
while keeping K = h cS k constant. 

In order to simulate a dissipativc environment in the 
quantum model we consider a master equation in the 
Lindblad form [28] for the density operator p of the sys- 
tem. We have 

2 2 

'p = -i[H s ,p] \ <T{Ll L„ p} + L.pLl (3) 
/i— l ii— l 

where H s = n 2 /2 + V(x,r) is the system Hamiltonian, 
are the Lindblad operators, and { , } denotes the an- 
ticommutator. We assume that dissipation is described 
by the lowering operators 

Li = gJ2 n VnTT\n) (n+l\, ( ) 

with n = 0, 1, ... [30]. Requiring that at short times (p) 
evolves like in the classical case, as it should be according 
to the Ehrenfest theorem, we obtain g = \J— In 7. 

The first two terms of Eq. (3) can be regarded as the 
evolution governed by an effective non-Hermitian Hamil- 
tonian, H cH = H s + iW, with W = -1/2 £ M L^Lfj,. In 
turn, the last term is responsible for the so-called quan- 
tum jumps. Taking an initial pure state \4>{tq)), the jump 
probabilities dp^ in an infinitesimal time dr are defined 
by dp^ = (</>(ro)|iJ l i M |^(ro))(ir, and the new states af- 
ter the jumps by |<^) = L^{t )) /\\L^{t ))\\. With 
probability dp^ a jump occurs and the system is left in 
the state With probability 1 — ^ dp^ there are no 
jumps and the system evolves according to the effective 
Hamiltonian H c g. In order to simulate the master equa- 
tion (3), we have used the quantum trajectories approach 
[31]. We start from a pure state |</>(to)) and, at intervals 
dr, we perform the following evaluation. We choose a 
random number e from a uniform distribution in the unit 
interval [0, 1]. If e < dp, where dp — tne s y s " 

tern jumps to one of the states |<^) (to |0i) if < e < dp\ 
or to |<^2) if dpi < e < dp\ + dpi). On the other hand, 
if e > dp, the evolution with the non-Hermitian Hamil- 
tonian H e g takes place, ending up in a state that we call 
\4>o). In both circumstances we renormalizc the state. We 
repeat this process as many times as n st0 ps = (Ar)/dr 
where At is the whole elapsed time during the evolu- 
tion. Note that we must take dr much smaller than the 
time scales relevant for the evolution of the open quan- 
tum system under investigation. In our simulations, dr 
is inversely proportional to 7j e ff- Given any observable 
A, we can write the mean value (A) T = Tr[Ap(r)] as the 
average over Af trajectories: 

<^>t=^^£<&0-)|4&(t)>. (5) 

i=l 



We found that a reasonable amount of trajectories (Af ~ 
100 — 500) is sufficient in order to obtain a satisfactory 
statistical convergence. We used up to 3 8 quantum levels 
to investigate the transition from quantum to classical 
behavior. 

The appearance of a strange attractor in our model is 
shown in the phase space portraits of Fig. 1, for param- 
eter values K = 7, 7 = 0.7, <j> = tt/2, a = 0.7 [32]. The 
left panel (classical mechanics) corresponds to a Poincare 
section constructed from 10 7 random initial conditions 
uniformly distributed in the area x <G [0, 2ir), p <G [— 7r; ir], 
after 100 iterations of map (2). From top to bottom we 
zoom in, showing the fractal structure of the attractor 
down to smaller and smaller scales. On the right column 
(quantum mechanics, ?i e fj = 0.012) we show the corre- 
sponding Husimi functions [33] averaged over iVj = 5 ini- 
tial conditions (randomly selected p values in the same 
interval as before), considering Af = 160 trajectories for 
each initial condition. We can see a good agreement be- 
tween the classical and the quantum phase space por- 
traits. Though the quantum version is less detailed the 
main classical patterns are well reproduced. On the other 
hand, a magnification of a small part of the phase space 
shows that the resolution of the quantum picture is lim- 
ited by the uncertainty principle. 

The attractor in Fig. 1 is strongly asymmetric, suggest- 
ing directed transport, that is, (p) ^ 0. This is confirmed 
by the numerical data of Fig. 2, where (p) is shown as a 
function of time, both in the classical and in the quantum 
case (for h cS = 0.037,0.11,0.33,0.99). In the classical 
case, (p) saturates after a very short time scale (approxi- 
mately 10 kicks) required for the setting in of the strange 
attractor. In the quantum case, a gradual approach to 
the classical limit can be clearly seen as h e g is reduced, in 
agreement with the correspondence principle The differ- 
ence S(p) between the classical and quantum saturation 
currents is shown function of fi c g in the inset of 

Fig. 2. In the same inset, we also show the numerical fit 
S(p) oc hps, with the fitting constant v 0.6. We see 
that the quantum fluctuations lead to a smaller value of 
the directed current as compared to the classical case. 

In Fig. 3, we show that it is possible to control the di- 
rection of transport by varying the phase <j). We observe 
that (p)-0 = —(p)<p- This is due to the fact that the po- 
tential V(x, t) has the symmetry V${x, t) — V-^—x, r). 
Therefore, the current can be reversed by changing <fi — > 
— 4>. In particular, there is a space symmetry for <fi = utt 
[in these cases, V(x,t) = V{— x,t)], which implies zero 
net current. This is confirmed, at </> = 0, from the nu- 
merical data of Fig. 3. 

We point out that quantum directed transport is ob- 
served in our model also for parameter values where the 
classical dynamics is characterized by a simple attractor 
(fixed points) instead of a strange one. Moreover, we have 
explored a different dissipation mechanism (i.e., loss of 
probability for \p\ larger than some threshold value pth 
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FIG. 1: (Color) Phase space pictures for K — 7, 7 = 0.7, 
4> — 7r/2, a — 0.7, after 100 kicks: classical Poincare sections 
(left) and quantum Husimi functions at fi c g = 0.012 (right). 
In the upper row, the displayed region is given by p G [—20, 20] 
and x £ [0, 27r) (note that, to draw the attractor, x is taken 
modulus 2tt). Magnifications of these plots are shown in the 
second and third rows (the area is reduced by a factor 1/9 and 
1/81, respectively). The color is proportional to the density: 
blue for zero and red for maximal density. 
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FIG. 2: Average momentum (p) as a function of time t (mea- 
sured in number of kicks), with same parameter values as in 
Fig. 1. The solid curve corresponds to the classical limit, 
while the other curves are quantum results at (from bottom 
to top) h e ff = 0.99,0.33,0.11,0.037. Each quantum curve is 
obtained from Ni = 60 initial conditions and M = 480 trajec- 
tories for each initial conditions. Inset: scaling of S{p) with 
h e ff (circles). The power- law fit 8{p) <x h" s gives v ~ 0.6 (full 
line) . 
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after each kick [34] ) and found that this case is also char- 
acterized by directed transport and a quantum strange 
repeller [35]. This is a further demonstration that quan- 
tum ratchets and asymmetric strange attractors and re- 
pellers are a typical feature of open quantum systems 
which are chaotic in the classical limit. 

It is important to analyze the stability of our ratchet 
phenomenon in respect to noise effects. For this pur- 
pose, we introduce memoryless fluctuations in the kick- 
ing strength: K — > K e (t) = K + e(t), where the noise 
values e(i) are uniformly and randomly distributed in 
the interval [— e, +e]. The dependence of the classical 



FIG. 3: Average momentum versus time, for (from top to 
bottom) cj> — 7r/2, 2n/5, 0, — 27r/5, — 7r/2, in the classical (full 
curves) and in the quantum case (dashed curves, R e ft = 0.11). 
The other parameter values are as in Figs. 1 and 2. 



and quantum currents on e is shown in Fig. 4. We can 
see that directed transport survives, both in the classi- 
cal and in the quantum case, up to a noise strength e of 
the order of the kicking strength K. We stress that the 
robustness of our ratchet model is in contrasts with the 
behavior of the quantum Hamiltonian ratchets discussed 
in Refs. [17, 18]. Indeed, noise eliminates the ratchet ef- 
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FIG. 4: Average momentum versus noise strength, in the clas- 
sical (full curve) and in the quantum case, for 7t e ff = 0.99 
(circles), 0.33 (squares), and 0.11 (triangles). The other pa- 
rameter values are as in Figs. 1,2. 



feet in Hamiltonian systems at large time scales, while 
in our case the ratchet is produced by dissipation and 
remains stable with respect to introduction of noise. We 
also note that the reduction of ratchet current with noise 
(Fig. 4) gives a qualitative reason for reduction of the cur- 
rent with increase of h e f / (Fig. 2) which leads to growth 
of quantum fluctuations. 

The model discussed here can be realized experimen- 
tally with cold atoms in a periodic standing wave of light, 
with parameters similar to those of Refs. [1, 2, 3, 4, 5]. 
For example, it is possible to use sodium atoms in a laser 
field with wave length Xl — 589 nm [1], which creates 
effective pulsed potential with classical chaos parame- 
ter K = \/2na$} c ffL} r T 2 and 7z e ff = 8uj r T. Here T is 
the pulse periodicity, Q e ff the effective Rabi frequency, 
uo r = hk\/2M the recoil frequency (fez, = 1/Az,, M 
atomic mass) and a is the fraction of Gaussian pulse 
duration in units of pulse period. For typical parameter 
conditions of Ref. [1], it is possible to take f2 e ff/27r = 75 
MHz, T = 0.8 /us, a = 0.05, which gives h cS w 1, K w 5, 
being similar to the case of Fig. 2 [36]. A dissipative fric- 
tion force can be created by additional laser fields which 
for sodium atoms can easily produce a dissipation rate 
2/3 w 4 x 10 5 s- 1 [37] giving 1 - 7 = 2/3T w 0.3. We also 
note that the Husimi function can in principle be mea- 
sured from a state reconstruction technique [38]. This 
could allow the experimental observation of a quantum 
strange ratchet attractor. 
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